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Abstract 

The quantum random walk has been much studied recently, largely due to its highly nonclassical 
behavior. In this paper, we study one possible route to classical behavior for the discrete quantum 
walk on the line: the presence of decoherence in the quantum "coin" which drives the walk. We 
find exact analytical expressions for the time dependence of the first two moments of position, and 
show that in the long-time limit the variance grows linearly with time, unlike the unitary walk. We 
compare this to the results of direct numerical simulation, and see how the form of the position 
distribution changes from the unitary to the usual classical result as we increase the strength of 
the decoherence. 
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I. INTRODUCTION 



In the classical discrete random walk, a particle is located at one of a set of definite 
positions (such as the set of integers on the line). In response to a random event — for 
example, the flipping of a coin — the particle moves either right or left. This process is 
iterated, and the motion of the particle is analyzed statistically. These systems provide 
good models for diffusion and other stochastic processes. 

Considerable work has been done recently on quantum random walks, which are unitary 
(and hence reversible) systems designed as analogues to the usual classical case. There are 
two general approaches to the problem: continuous [3, 0, HI and discrete [H, il, El 0, IE 
HQBHQEHEHBIIIHB unitary walks. This paper is exclusively 
concerned with the discrete walk. In this discrete case, we introduce an extra "coin" degree 
of freedom (usually a single quantum bit) into the system. Just as in the classical random 
walk, the outcome of a "coin flip" determines which way the particle moves; but in the 
quantum case, both the "flip" of the coin and the conditional motion of the particle are 
unitary transformations. Different possible classical paths can interfere with each other. 

In this paper we look at quantum walks on the infinite line. The particle is initially at 
position x = and is free to travel off to infinity in either direction. We will look at both the 
probability distribution in p(x,t) = (x\pt\x), and at the long-time behavior of the moments 
(x) and (x 2 ) — (x) 2 as functions of t. 

For a classical random walk, p(x, t) has the form of a binomial distribution, with a width 
which spreads like the variance x 2 — x 2 grows linearly with time. The variance in the 
quantum walk, by contrast, grows quadratically with time; and the distribution p(x,t) has a 
complicated, oscillatory form. Both of these are effects of interference between the possible 
paths of the particle. 

It should be possible to recover the classical behavior as some kind of limit of the quantum 
system. There are two obvious ways to regain the classical result. If the quantum "coin" is 
measured at every step, then the record of the measurement outcomes singles out a particular 
classical path. By averaging over all possible measurement records, one recovers the usual 
classical behavior [ill ]. 

Alternatively, rather than re-using the same coin every time, one could replace it with a 
new quantum coin for each flip. After a time t one would have accumulated t coins, all of 
them entangled with the position of the particle. By measuring them, one could reconstruct 
an unique classical path; averaging over the outcomes would once again produce the classical 
result. 

These two approaches, which are equivalent in the classical limit, give two different routes 
from quantum to classical [23] . We might increase the number of coins used to generate the 
walk, cycling among M different coins, in the limit using a new coin at each step. Or 
we might weakly measure the coin after each step, reaching the classical limit with strong, 
projective measurements. This is equivalent to having a coin which is subject to decoherence. 



In another paper |2^ we have considered the quantum random walk with multiple coins. 
In this case, the quantum behavior remains qualitatively unchanged until we reach the limit 
of a new coin for each step, at which point classical behavior is recovered. 

In this paper, we consider the quantum random walk with a single coin subject to deco- 
herence. We will see that in this case, the behavior is qualitatively quite different from the 
unitary quantum random walk. The usual classical solution is recovered in the limit where 
the coin decoheres completely every step; but even with weaker decoherence, the variance 
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of the position distribution grows linearly with time, rather than quadratically (as in the 
unitary case). 

In section II we present an analytical result for the moments of the decoherent walk, and 
compare them to the results from direct numerical simulations. In section III, we see how 
the probability distribution p(x, t) changes as we introduce decoherence. Finally, in section 
IV we summarize our results and state conclusions. 

II. MOMENTS OF THE DECOHERENT WALK 
A. The unitary walk on the line 

Let us now consider a fairly general quantum random walk on the line. The particle 
degree of freedom has a basis of position eigenstates {\x)} where x can be any integer. The 
position operator is x, and x\x) = x\x). We will assume that the particle begins the walk 
at the origin, in state |0). The walk is driven by a separate "coin" degree of freedom: a 
D-dimensional system with an initial state |$ ). Let Vr : Vi, be two orthogonal projectors 
on the Hilbert space of the "coin," such that + Vl — I- These represent the two 
possible outcomes of the coin flip: Heads or Tails, Right or Left. We also define a unitary 
transformation U which "flips" the coin by rotating a coin showing heads or tails into a 
superposition of the two. One step of the quantum random walk is given by the unitary 
operator 




(HI) 



where S, are unitary shift operators on the particle position: 




(II.2) 



The full initial state of the system (particle and "coin" ) is 




(II.3) 



We can identify the eigenvectors 





(II.4) 



X 



with eigenvalues 



S\k) 
&\k) 




(II.5) 



The inverse transformation is 




(II.6) 



In particular, the initial state of the particle is 




(II.7) 
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These state vectors \k) are not renormalizable, but if used with caution they greatly simplify 
the calculations. In the k basis, the evolution operator becomes 



E (\k) <g> |$)) = \k) ® (e- lk V R + e< fe P L ) U\$) , 

= \k)®U k \$) , (II.8) 

where Uk is also a unitary operator. 

The usual case considered in the literature has taken the coin to be a simple two-level 
system, and the "flip" operator U to be the usual Hadamard transformation H: 

H\R) = -L(|i2) + |L» , 

H\L) = -L(\R)-\L}) . (II.9) 

The projectors are "Pr = \R)(R\, Vl = \L)(L\. The walk on the line in this case has been 
exactly solved by Nayak and Vishwanath p. 

For the present, we will continue without assuming a particular form for U, Vr or PL- 
Later we will specialize to make comparison to numerical simulations. 



B. Decoherence 

We now generalize to allow for decoherence. Suppose that before each unitary "flip" of 
the coin, a completely positive map is performed on the coin (note, NOT on both the coin 
and the particle). This map is given by a set of operators {A n } on the coin degree of freedom 
which satisfy 

Y,AiA n = I. (11.10) 

n 

A density operator x f° r the coin degree of freedom is transformed 

X^X' = J2 A ^ A n- ( IL11 ) 

n 

If we apply this to a general density operator for the joint particle/coin system 

P = /S/ f^^'i®^'' ( IL12 ) 
then after one step the state becomes 

Y \ ^\k){k'\®Y,U k A nXkk 'Alul . (11.13) 

The initial state is 

po = /§/ d ^\k){k'\®\^)(M- (n.i4) 
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Let the quantum random walk proceed for t steps. Then the state evolves to 

/dk f dk 1 
^ —\k)(k>\® ^•■•^n 1 |$0)($0|^ 1 ^V--il^- (H.15) 

ni,...,n t 

We define the superoperator C kk i on the coin degree of freedom 

/§/ ® Y,U k A nXkki mi = J^J ^\k){k'\®£k k 'Xk« ■ (H.16) 

In terms of the superoperator, 

/dk C dk' 
f£ J ^\k)(k'\®jCU$ )(%\. (IL17) 

Note that for k = k' this superoperator preserves the trace. This implies that 

Tr(££ fc d} = Tr{6) , (11.18) 



for any operator O. This identity will prove useful later. 
The probability to reach a point x at time t is 

p(x,t) = Tr {\x) (x\p t } = (x\p t \x) 
1 



dk / dk'{k\x)(x\k')Ti {Ci k/ \%}(%\} 



(27T) 2 

= j^rl dk l ^ ,e "" (fc " fc ' )Ti '{4 fc ^o)($oi} . (ii.i9) 

C. Moments of position 

Eq. (jII.19|) for p(x, t) will be difficult to evaluate, in general. However, we can get con- 
siderably further by restricting our interest to the moments of this distribution. 

(x m ) t = 5> to pGM) 

X 

= J^fY.*™ j dk j dk'e-^ k - k '^r{Cl k ,\%}(%\} . (11.20) 

We can then invert the order of operations and do the x sum first. This sum can be exactly 
carried out in terms of derivatives of the delta function: 

— x m e- ix{k ~ k>) = (-i) m 5 {m \k - k') . (11.21) 

X 

Inserting this result back into our expression for (x m ) t yields 

(x m )t = j dkj dk'^ m \k - fc')Tr {4 fe ,|$ )($o|} . (11.22) 



We can then integrate this by parts. 

In carrying out this integration by parts, we will need 



d 

dk 



Tr 



{Afe'O} = Tr {(dC kk ,/dk)6 

= J2 TT {( d Uk/dk)A n dAlul 



[11.23) 



where 



d£k 
dk 

dUl 
dk 

z 



-i{Vo-V 1 )U k = -iZU k 



iulz 



v R -v L 

Substituting this back into (|II.23J) we get 

d 



(11.24) 



dk TT { Ckk '° 



-iTi{{C kkl 6)Z 

-^ Tr { Ckk ' d 



(11.25) 



Making use of (|II.18j) and (|II.25jl . when we carry out the integration by parts for the first 
moment we get 



3=1 J 



We can simplify our notation slightly by defining 

We can carry out a similar integration by parts to get the second moment: 

t j 



(11.26) 



(11.27) 



. dk 

2tt 

t 3-1 



.3=1 j'=l 



+ 



^^^{z£r'((4i$ )($oi)^) 



3=1 j'=l 



(11.28) 



Note that this form is rather similar to that for a correlation function. 

For the unitary walk (23l . 24]|, we were able to separate the expressions for the moments 
into oscillatory and nonoscillatory terms, and consider only the nonoscillatory terms for 
the long-time limit. In the present case, this will not work; because of the presence of 
decoherence, the oscillations are damped. However, the long-time limit still simplifies in 
another way. If we think of C k as a linear transformation, all its eigenvalues must clearly 
obey |A| < 1 in order for it to be completely positive. In finding the long-time limit, we 
need pick out only those components of the expressions (jll.26j) and (jll.28|) which do not die 
away at large t. 
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D. The Hadamard walk with decoherence 



To get any further than this, we need to specialize to a particular model. Let's choose the 
standard two-dimensional "coin" and the Hadamard walk as described above. This makes 
our operator U k 

ft = ^ ( e « _ e «) • (H.29) 

We need also to pick a particular form for the decoherence. I will choose the decoherence 
produced by the three operators 

A = y /p\R)(R\ , 
A 1 = ^\L)(L\ , 

M = ^fl-pl ■ (H.30) 

This resembles a coin which has a probability p per step of being measured; however, this 
is equivalent to a large class of other decoherence models, such as pure dephasing 

A' = ^=(e^\R)(R\+e^\L)(L\) , 

A l = ^{e^\R){R\+e^\L){L\) , (11.31) 

or weak measurement 



AZ = y/q\R)(R\ + y/l=l\L)(L\, 

A'{ = y/T^\R){R\ + y/q\L)(L\ . (11.32) 

All three of these decoherence processes represent the same completely positive map on 
the density matrix; one can transform from one to the other by invoking the relationships 
between the parameters 

1 - p = cos 2^ = 2y/q{l-q) . (11.33) 

We will generally use the first in our derivations, for simplicity; but the second is more 
convenient for numerical calculations. Note that this equivalence applies to the average 
behavior over all outcomes; for a single sequence of outcomes, the conditional dynamics of 
these three decoherence (or weak measurement) models can look quite different. 

In the absence of decoherence, the Hadamard walk exhibits a linear drift in (x) and a 
quadratic growth in the variance (x 2 ) — (x) 2 . We plot these moments in figure 1. 

Because C k is linear, we can represent it as a matrix acting on the space of two-by-two 
operators. A convenient representation is to write 

O = r I + r\0\ + r 2 cr 2 + r 3 a 3 , (11.34) 

where Ci,2,3 = a x,y,z are the usual Pauli matrices. We can then represent O by a column 
vector 



o 



7-1 



(11.35) 
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(a) Co) 

FIG. 1: The (a) expectation (x), and (b) variance (x 2 ) — (x) 2 , for the Hadamard walk on the line 
without decoherence. The coin begins in state \R). Note that the first moment exhibits a linear 
drift with time, which seems to reflect a peculiar "memory" for the initial state; a system starting 
with the coin in state \L) would drift symmetrically in the other direction. The variance grows 
quadratically in time, in contrast to the linear growth in the classical random walk. Note that if 
viewed at a finer scale, both of these curves would exhibit oscillations about the simple power-law 
behavior. 



The action of Ck on O is then given by 



c k o 



A o 



\0 l-p 



\ 

[1 — p) sin2k cos2fc 
'1 — p)cos2k — sin2/c 
J 



M 

ri 

w 



(11.36) 



Note that since r = Tr{0}, it is unaffected by Ck, which is trace-preserving. So the only 
nontrivial dynamics result from the three-by-three submatrix 



M, 




(1 — p) sin2fc cos2/c 
1 — p) cos 2k —sin 2 A; 




(11.37) 



We also need to know the effects of left-multiplying and right-multiplying by Z. These 
are given by the two matrices 



Zr = 



(0 








1\ 








i 








—i 








V 











Z R = 



/o l\ 
-i 
i 

\1 0/ 



(11.38) 



Finally, taking the trace picks out the component of the column vector and drops the rest. 

Let us take these expressions and apply them to equation (jll.26|) for the first moment. 
In the integrand, the initial density matrix for the coin is multiplied j times by Ck, then 
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left-multiplied by Z, and finally the trace is taken. Given the above expression for Z^, we 
see that this is the same as multiplying the three- vector (ri,r2,r^) j times by M k and then 
keeping only the r% component of the result. This gives us the new expression: 



(x) t = -1. j dk (0 1) 



U=i 




= -^f dk(0 1) [(l-M k )-\M k -M t k +1 )} ^r 2 j . (11.39) 

The eigenvalues of M k are complicated, but fortunately we don't need to evaluate them. All 
we need to know is that all of them obey < |A| < 1 (where both these inequalities are 
strict). In the long time limit, therefore, M l k +1 — > 0, and the moment becomes approximately 

(x) t ^-^Jdk(00 1) [(1 - M k )- l M k ] \ rz j . (11.40) 

Note that all t dependence has vanished. Therefore, in the long time limit, the first moment 
tends to a constant. 

As it happens, the matrix 1 — M k is exactly invertible: 

j / 1 + (1 — p) cos2A; — (1 — p)sin2fc (1 — p) + cos2/c 

(\-M k y 1 = — -(l-p)sin2Jfe 1- (l-p)cos2Jfe -sin2A; 

P( 2 -P) \(1 - p )(l + (1 -p)cos2k) -{l-p) 2 sin2k 1 + (1 - p) cos 2k / 

(11.41) 

Inserting this and (jll.37|) into (jll.40|) yields 
If 1 -p 

(x) t w — — J dk~— — —- [(1 - p)r 3 + ri + (r 3 + (1 — p)r 1 ) cos 2k — (1 — p)r 2 sin 2k] 
1 — p 

[(1 -p)r 3 + n) 



p{2-p) 

1 — p 
p(2-p) 



[(l-p)(\a\ 2 -\(3\ 2 ) + (a*(3 + a(3*)} , (H.42) 



where the initial state of the coin is |$o) = ot\R) + (3\L). 

In figure 2, we compare this result to the results of direct numerical simulation for the 
initial state |$o) = \R)- As we can see, the first moment does tend to drift towards a 
constant value asymptotically. What variation there is is a result of statistical error in the 
Monte Carlo simulation rather than a poor fit. 

The second moment is a somewhat more complicated calculation. Let us rewrite the 
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6 



5 
4 



a 3 

v 



Analytical 

6=71/16 
6=71/8 
6=37t/16 
6=7c/4 






50 100 150 200 250 300 350 400 450 500 

t 

FIG. 2: (x) vs. t for the Hadamard walk on the line with decoherence, for 9 = 
7r/16, it/8, 37r/16, 7t/4. For all cases the coin began in the initial state \R). Note that (x) goes 
asymptotically to a constant value at long times which matches our analytical estimate; this drift 
goes to zero with increasing decoherence, vanishing at 8 = ir/4 (i.e., p = 1). Note that the 
irregularities in the broken curves reflect statistical errors in the Monte Carlo calculations. 



equation (jll.28|) in terms of our four-by-four matrices: 

1 r " EE IY {^r i '(^4>o)($oi) 
u=i j'=i 

+EE IY { i£ r i '((4>o)(*oi)^) 

j= i j'=i 



(x ), = / dk 



t - — I dk (1 0) 



/ - — / dk (0 1) 



/1\ 

n 

\r 3 J 



.i=i j'=i 



/1\ 

\r 3 / 



(11.43) 



The initial t term comes from the j = j' components of (jll.28|) . where the two factors of Z 
cancel out. 
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Recall the block- diagonal form (|II.36|) of Ck, and note that 



{Z L + Z 



R 



/0 2\ 





\2 o o oy 



(11.44) 



is extremely sparse. In this case, it makes sense to separate the and 1,2,3 components of 
the expression: 



(x 



j=l j'=l J 



cr (z L + Zr )c{ 



/A 

o 
o 

w 



+c{- 3 \z L + z R )c 



n 

w 



(11.45) 



We can drop the second term, because it will be cancelled in the inner product with 
(0 l). This means that (x 2 ) t has no dependence on the initial state! The remaining 
term becomes 



*-££^/^( 1 )£r' 

7 = 1 j' = l J 



t - 



t - 



1 

2^ 
1 

27 



d£; (0 1) 
dk (0 1) 



t j-i 



.i=i i'=i 



M 

o 



W 



x [t - (1 - M k )~ x M k + (1 - Mfc)- 1 !/^] (1 - M fe ) _1 M fe ( 



(11.46) 



In the long-time limit, we can drop the M l k terms as negligible. We can thereby substitute 
the exact expressions (jll.37|) and (jll.41|) for the remaining matrices and simplify: 



(x 2 ). — t H / dk— — — - 

V h 2nJ p{2-p) 



2i(l -p + cos2fc) 



1 



p(2-p) 

2(1 -P) 2 



t[l + 



(4(2 - 2p + p 2 ) cos 2k + (1 - p) (7 + cos 4Jfe)) 
7(1 -P? 



p{2 — p) J p 2 (2 — p) 2 



(11.47) 
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Analytical 

6=7t/16 

e=7t/8 
e=37r/i6 
e=ic/4 



20 



-4- 



2000 
1000 


50 100 150 200 250 300 350 400 450 500 

t 

FIG. 3: (x 2 ) — (x) 2 vs. t for the quantum random walk with decoherence, for 6 = 
vr/16, vr/8, 37r/16, 7r/4. For all cases the coin began in the initial state Note that the vari- 
ance goes asymptotically to a linear growth at long times which matches our analytical estimate; 
the rate of growth goes to one with increasing decoherence, matching the classical case at 6 = 7r/4 
(i.e., p = 1). The inset shows the time-dependence of the variance at short times, where it still 
exhibits quadratic growth. 



Thus, in the case of the decoherent coin the variance will grow linearly with time at long 
times, just as in the classical case, though the rate of growth will be greater than for the 
classical random walk. Note that as p — *■ 1, (x 2 ) t — > t, which is the classical result. 

Because of the dependence of the first moment on the initial condition, the variance will 
have some residual dependence on the initial condition; but since this is only a constant it 
is quite unimportant. Upon comparison to the numerical results (see figure 3), we see that 
this analytical expression for the variance matches extremely well. 

Since this approximation holds in the long-time limit, it is reasonable to ask how long is a 
"long time." For weak decoherence (p <C 1), the evolution exhibits the same behavior as the 
unitary walk for short times, with the second moment growing quadratically in time, and 
only switches over to linear growth past some critical timescale. This changeover is visible 
in figure 3. 

We take the long-time approximation by neglecting the M l k term in (jll.46j) . If we examine 
M|, we see that it can be written as (1 — p) multiplying a matrix whose determinant is less 
than one. Therefore, if we want this term to be negligible, then 

\\Ml\\ < (1 -p) t/2 < e (11.48) 
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for some positive e C 1. This implies that 

p > (2/t) log(l/e) = c/t (11.49) 

which implies in turn that the second moment of position will grow linearly at t > c/p, for 
a large enough constant c. For t < c/p there will be a transition over to quadratic time 
dependence. For small p, equation (jll.47|) becomes 

(x 2 ) t ==t(l + l/p) - 7/Ap 2 . (11.50) 

III. THE POSITION DISTRIBUTION 

As mentioned before, Eq. (jII.19|) for p(x, t) is difficult to evaluate analytically. It is 
straightforward, however, to solve this system numerically. Rather than solve for the den- 
sity matrix p, we instead have used quantum trajectory techniques |25[ to do a quantum 
Monte Carlo simulation, averaging over many runs to find both the distribution itself and 
its moments. 

For the purposes of these simulations, it proved more convenient to use the form (|11.31|) for 
the decoherence process. We have therefore labeled the figures with the appropriate values 
of the dephasing parameter 9 instead of p. In performing the simulations, after each step of 
the quantum walk we randomly applied either operator a/2Aq or V2A[ from eq. (jll.31|) to 
the state, choosing them with equal probability. Because V%A' 1 are unitary operators, it 
was unnecessary to renormalize the state. We then found the mean probability distributions 
p(x) by averaging over many runs; each distribution plotted in this paper represents an 
average over 10000 runs of the Monte Carlo program. In figure 4 we see how the probability 
distribution changes as we go to the classical limit of complete decoherence at every step 

(0=7r/4). 

As we see from the figure, the presence of decoherence quickly wipes out the most con- 
spicuous signs of interference: the peaks at ±t/^/2 and the oscillations in the distribution. 
A new peak, centered at t — 0, makes its appearance, though for weak decoherence it does 
not have the usual binomial form. However, some residual effects of interference persist up 
until just short of the classical limit. In particular, the distribution is broadened, with long 
"tails," and displaced somewhat from the central position. This is reflected in our long-time 
solution for the variance (|II.47|) . which is higher in the quantum case than the classical 
except in the limit p — > 1. 

In the unitary case, interference effects cause the evolution to "remember" the starting 
state in a non-intuitive way. As we saw in figure 1, if we start the coin in state \R), the 
system retains a linear drift to the right for all time. If we had instead started in the state 
\L), the system would have drifted symmetrically to the left. This "memory" also makes 
the position distribution asymmetric, as we see in fig. 4a. 

The presence of decoherence eliminates this effect at long times: the system "forgets" 
the initial conditions after a while. This effect will still persist at short times, however, and 
produce a tendency for the particle to move repeatedly in the same direction. No doubt this 
tendency produces the broadening in the position distribution and the significant tails. 
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IV. CONCLUSIONS 



We have examined a possible route from quantum to classical for the quantum random 
walk: letting the quantum coin undergo decoherence with time. Using the long-time behavior 
of the position moments as a qualitative marker of classicality, we see that even very weak 
decoherence changes the growth of the variance from a quadratic to a linear function of time. 
Since the usual classical random walk has a linear growth of the variance, in a particular 
sense we can claim that the decoherent walk is indeed classical. 

This situation is quite different from the quantum random walk with multiple coins 24]| ; 
for that system, quadratic growth remained except in the limit of a new coin for every step. 
One might reasonably claim that the multicoin system remains "quantum" even in the limit 
of very large numbers of coins, while the decoherent system remains "classical" even in the 
limit of very weak noise. 

In spite of this, some effects of interference are important even in the presence of deco- 
herence. In particular, the variance grows more rapidly in the quantum than the classical 
case. This may reflect a tendency for the particle to be "biased" in one direction for short 
times, while interference effects enable it to "remember" its starting state. This is consistent 
with the position distributions, which are broader than the classical binomial distribution, 
and have nontrivial tails. 

In the long time limit, the time dependence of the moments becomes tractable because 
the evolution superoperator has only one eigenvalue of modulus 1. All other components 
decay away exponentially with time. It is interesting to speculate whether a decoherence 
process with several modulus 1 eigenvalues might continue to exhibit quantum behavior. 
This is probably possible, at least for higher-dimensional coins. 

We should emphasize that the system we have studied in this paper has decoherence only 
in the coin degree of freedom; any effects of decoherence on the particle are indirect. It might 
be arguably natural to include decoherence in the particle position as well. Such systems 
have been studied numerically, and exhibit interesting effects of their own j^"" 
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FIG. 4: The probability distributions p(x, t) at t = 500 for the quantum random walk with 
decoherence: (a) 9 = (unitary case), (b) 9 = tt/16, (c) 9 = tt/8 (d) 9 = 3vr/16, (e) 9 = vr/4 
(classical case). A central peak appears when decoherence is included, which becomes increasingly 
dominant; however, the peak is broadened compared to the classical case, with nontrivial tails 
which disappear only in the classical limit, when p(x, t) goes over to the binomial distribution. 
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